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Abstract. Dissipation and fluctuations of one-body observables in heavy-ion reactions around the Coulomb 
barrier are investigated with a microscopic stochastic mean-field approach. By projecting the stochastic mean- 
field dynamics on a suitable collective path, transport coefficients associated with the relative distance between 
colliding nuclei and a fragment mass are extracted. Although microscopic mean-field approach is know to un- 
derestimate the variance of fragment mass distribution, the description of the variance is much improved by the 
stochastic mean-field method. While fluctuations are consistent with the empirical (semiclassical) analysis of the 
experimental data, concerning mean values of macroscopic variables the semiclassical description breaks down 
below the Coulomb barrier. 



1 Introduction 

The interplay between nuclear structure and dynamical ef- 
fects is crucial to properly describing heavy-ion fusion re- 
actions at energies close to the Coulomb baiTier The mean- 
field theory based on the Skyrme energy density functional 
provides a rather unique tool for describing nuclear struc- 
ture and nuclear reactions over the whole nuclear chart in 
a unified framework. For nuclear dynamical calculations, 
time-dependent Hartree-Fock (TDHF) model has been de- 
veloped [1 2 3 4|. This model automatically includes im- 
portant dynamical effects such as vibration of nuclei, neck 
formation, and nucleon transfer during reactions. Most re- 
cent TDHF simulations are able to include all terms of the 
Skyrme energy density functional, for example the spin- 
orbit terms, and to break spacial symmetries |5 6 7 8 9 101. 
To perform full three-dimensional calculations and to use 
effective forces consistent with static calculations is crucial 
to accounting for the richness of nuclear shapes in dynam- 
ical evolution. 

In the mean-field theory, short range two-body correla- 
tions are neglected and nucleons move in the self-consistent 
potential produced by all other nucleons. This is a good 
approximation at low energies since Pauli blocking is ef- 
fective for scattering into unoccupied states. Consequently, 
collective energy is converted into intrinsic degrees of free- 
dom via interaction of nucleons with the self-consistent 
mean field, so-called one-body dissipation. One-body dis- 
sipation mechanism plays a dominant role in low energy 
nuclear dynamics. One important limitation of the mean- 
field theory is related with dynamical fluctuations of col- 
lective variables. In the mean-field description, while single- 
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particle motion is treated in a quantal framework, collec- 
tive motion is treated almost in a classical approximation. 
Therefore, TDHF provides a good description for average 
evolution; however it severely underestimates fluctuations 
of collective variables 01 II12I13I . 

After the first application of the TDHF theory, much 
eff'ort has been devoted to overcoming this difficulty and to 
developing transport theories that are able to describe not 
only mean values but also fluctuations (for a review, see 
Refs II14I15I ). Among them, the variational principle by 
Balian and Veneroni appears as one of the most promis- 
ing methods 1 16il7J . However, even nowadays it remains 
difficult to apply II18II9I . The absence of a practical so- 
lution to include fluctuations beyond mean field in a fully 
microscopic framework strongly restricts applications of 
mean-field-based theories to low energy nuclear reactions. 

A stochastic mean-field (SMF) approach has been pro- 
posed for describing fluctuation dynamics |20|. For small 
amplitude fluctuations, this model gives a result for disper- 
sion of a one-body observable that is identical to the result 
obtained through a variational approach 1161 . 

In this contribution, first, we discuss the property of 
nucleus-nucleus potential and one-body dissipation deduced 
from TDHF calculations. They are not constrained by adi- 
abatic or diabatic approximation. Therefore, this method 
should provide an accurate description of nucleus-nucleus 
potential and one-body dissipation at energies around the 
Coulomb barrier [21 22 1. Then, we discuss how to over- 
come the failure of description of fluctuations of one-body 
observables in TDHF and employ the SMF approach II20I23I 
We project the SMF evolution on a collective path to ob- 
tain the expression of diffusion coefficients. We show that 
the description of the variance of fragment mass distribu- 
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tions in transfer reactions by the stochastic mean-field ap- 
proach is much improved compared to that of the TDHF 
model l,24J. 



2 Nucleus-nucleus potential and one-body 
dissipation from mean-field dynamics 

Nucleus-nucleus potential and one-body dissipation are ex- 
tracted as follows |21 22 1 : (i) We solve the TDHF equation 
for head-on collisions to obtain the time evolution of the 
total density of colliding nuclei, (ii) Defining a separation 
plane between two nuclei, we compute at each time the 
relative distance R, associated momentum P, and reduced 
mass fi of colliding nuclei, (iii) We assume that mean-field 
evolution obeys a classical equation of motion with a fric- 
tion term: 
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Fig. 1. Barrier height Vg deduced from TDHF minus experimen- 
tal barrier height Vg"'' as a function of extracted barrier height for 
the reactions indicated in the figure. The values Vb are deduced 
from high energy (solid line) and from low energy (Ecti. ~ Vg) 
(dashed line) TDHF trajectories, respectively. 
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where V(R) and denote the nucleus-nucleus poten- 
tial and friction coefficient, respectively. The friction coef- 
ficient yiR) describes the effect of one-body energy dissi- 
pation from the macroscopic collective degrees of freedom 
to the microscopic ones. For the TDHF calculations pre- 
sented in this contribution, the three-dimensional TDHF 
code developed by P. Bonche and coworkers with the SLy4d 
Skyrme effective force [5 1 is used. The mesh sizes in space 
and in time are 0.8 fm and 0.45 fm/c, respectively. For 
more details, see Refs. 02 11221 . 

Figure[T]shows the difference between the barrier height 
deduced from TDHF evolution {Vb) and the experimental 
barrier height (V^"'') 125] as a function of Vb- The solid 
line corresponds to the barrier height extracted using high- 
energy TDHF trajectories (Zscm. » Vb), whereas the dashed 
line is the result for low-energy TDHF trajectories (Ec.m. ~ 
Vb)- The former identifies with the barrier height of the 
frozen density approximation |26|. Dynamical reduction 
of the barrier height from high-energy TDHF to low-energy 
TDHF is clearly seen for all reactions. Moreover, because 
of this reduction, the value of the barrier height at low en- 
ergy approaches the experimental data. This underlines the 
importance of dynamical effects close to the Coulomb bar- 
rier and shows the precision of our method. 

This method also provides information on one-body 
dissipation through the friction coefficient y. In Fig.|2] we 
present reduced friction coefficients /3{R) = y{R)lii as a 
function of R scaled by the Coulomb barrier radius Rb for 
different systems. This clearly shows that the order of mag- 
nitude of [i(R) and the radial dependence are almost inde- 
pendent on the size and asymmetry of the system. We com- 
pare these friction coefficients with that of a microscopic 
model based on small amplitude response by Adamian et 
al. [27 1 by the solid circles. The radial dependence and the 
magnitude of their friction coefficient are very similar to 
those extracted from TDHF evolution. 
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Fig. 2. Extracted reduced friction p{R) = y(R)/fi as a function 
of R scaled by the Coulomb barrier radius Rb for difi'erent sys- 
tems. A microscopic friction coefficient by Adamian et al. 1 27 J is 
shown by the solid circles for comparison. 



3 Fluctuations of one-body observables 
from stochastic mean-field dynamics 

In the previous section, we show that mean-field dynamics 
gives a good description for the nucleus-nucleus potential 
and one-body dissipation. However, it is well known that 
TDHF can not reproduce fluctuations of one-body observ- 
ables, for example the variance of fragment mass distribu- 
tions in deep inelastic collisions lfT3l . although TDHF cal- 
culations well reproduce the mean value of fragment mass. 

In order to overcome this difficulty, we have recently 
proposed a stochastic mean-field (SMF) approach, which 
is a stochastic extension of the mean-field model for low 
energy nuclear dynamics so as to include zero-point fluc- 
tuations of the initial state |20 23 1. The initial density fluc- 
tuations are simulated by representing the initial state in 
terms of a suitable ensemble of initial single-particle den- 
sity matrices, which is similar to the idea in Refs. II28I29I . 
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In this manner, the description with a single Slater deter- 
minant is replaced by a superposition of several Slater de- 
terminants. A member of the ensemble of density matrices, 
indicated by event label A, can be expressed as 



,(r,t;A)pUcrT)0j^r(r',t;A), (2) 



IJtTT 



where the sums / and j run over a complete set of single- 
particle wave functions ?&,o-T(r, t; A) with spin-isospin quan- 
tum numbers cr, r. According to the description of the SMF 
approach Il20l . matrix elements p'^-icn) are assumed to be 
time-independent random Gaussian numbers with a mean 
value pfjicTT) = Sijti""^ and with a variance 



6pfj(cTT)6pj.,(cr'T') 
1 



where X denotes the ensemble average of X. Here, n""^ de- 
notes the average single-particle occupation factor At zero 
temperature occupation factors are and 1, and at finite 
temperature they are determined by the Fermi-Dirac distri- 
bution. The great advantage of the SMF approach is that 
each Slater determinant evolves independently from each 
other following the time evolution of its single-particle wave 
functions in its self-consistent mean-field Hamiltonian, de- 
noted by hip'^), according to 

m-t>i^riX, t- A) = hifi^)t>i^r(r, t; A). (4) 
at 

In the following applications, we focus on the head- 
on collision along the x axis around the Coulomb barrier 
energy. 



3.1 Fusion reactions 

First, we apply the SMF approach to fusion reactions | 
We project the SMF evolution on a one-dimensional macro- 
scopic Langevin equation, which is similar to Eq. ([TJ ex- 
cept an additional Gaussian random force ^p{t): 



dR 



(5) 



Ignoring non-Markovian effects, the random force ^p(f) 
with zero mean value reduces to white noise specified by 
the following correlation function: 



epmiit') = 2s{t - t')Dpp(R). 



(6) 



Here Dpp(R) denotes the momentum diffusion coefficient. 
Denoting xq as the position of the separation plane between 
the two nuclei (xq = in this case), we have the following 
semiclassical expression for the nucleon diffusion coeffi- 
cient according to Ref . Il23l : 
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Fig. 3. Diffusion coefficient as a function of the reiative distance 
for the head-on '*''Ca+'*"Ca coliision at „ = 100 MeV. 



J 2nti m 2 

X 2{/r(^o, P., f) [1 - /r(^o, p., t)l Q{xa, t)] 

(TT 

+ f^i^o. P., [1 - /r(^o, p., t)] /Q(xo, t)}. (7) 

Here m is the nucleon mass, Q(xo, f) is ffie phase-space vol- 
ume over the window given in Ref. 1231 . and 



fpJj{XQ,Px, t) 

■■ JJ dydz J ds^exp^-py.Sj.'j 



ieP/T 



(8) 



is the average value of ffie reduced Wigner distribution 
associated with single-particle wave functions originating 
from the projectile/target. We note that the expression of 
the diffusion coefficient (|7]) has the same form as that given 
by the phenomenological nucleon exchange model in Ref. 
Il30l . We also note that diffusion coefficients can be evalu- 
ated in terms of the average TDHF evolution through the 
Wigner transformation dS). 

As an example, the momentum diffusion coefficient for 
the head-on ''"CaH-''"Ca collision at fi^.m. = 100 MeV is 
shown in Fig. [3] 



3.2 Variance of fragment mass distribution in 
transfer reactions 

We show another application of SMF to the variance of 
fragment mass distribution in transfer reactions to improve 
the mean-field description Il24l . To do so, we investigate 
head-on transfer reactions at energies just below the Coulomb 
barrier, where nucleon exchange will occur during reac- 
tion, and estimate the variance of fragment mass distribu- 
tion. 

By projecting the SMF evolution on the collective space 
for the fragment mass number, time evolution of the mass 
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number of the projectile-like fragment Ap for an event A is 
described by a Langevin equation ||3T1 . 

jA'p^u(Aj„t)+^'^(t), (9) 

where v(Ap, f) denotes the drift term for nucleon transfer. 
The Gaussian white noise random force ^'\{t) is determined 
with zero mean value and a correlation function. 



f;J(f)^j(f') = 25(f-f')£'AA, (10) 

where Daa is the diffusion coefficient associated with nu- 
cleon exchange. The variance cr^^ of fragment mass dis- 
tribution is determined by small fluctuations of the mass 

number 5Ap through cr^^if) - SApSAp. The expression of 
the diffusion coefficient Daa is similar to the momentum 
diffusion coefficient Dpp except that p^. is removed from 
the integral in Eq. (|7]i. Again, we note that the nucleon dif- 
fusion coefficient can also be evaluated from the average 
TDHF evolution. 

According to the Langevin equation, neglecting contri- 
butions from the drift term, the variance is related to the 
diff'usion coefficient according to 1311321 

crii(f)-2 r DAA(s)ds. (11) 
Jo 

In the phenomenological nucleon exchange model, the re- 
lation cr^yi(r) = A^exc(0 was obtained, where A^exc(0 denotes 
the accumulated total number of exchanged nucleons until 
time t, and was extensively used to analyze the experimen- 
tal data of mass variance |33 34|. In the following, to check 
whether the SMF approach satisfies the above relation, we 
estimate the both quantities by the SMF approach. 

We carry out calculations for the head-on "^"Ca-i-^^Ca 
reaction at energies just below the Coulomb barrier Fig- 
ure ID illustrates the dependence of diffusion coefficients 
at different center-of-mass energies. The Coulomb barrier 
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Fig. 4. Time evolution of diffusion coefficient calculated in the 
SMF approach for the head-on "^''Ca + "^''Ca reaction at different 
center-of-mass energies below the Coulomb barrier. 
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Fig. 5. Time evolution of cr^^ obtained from the SMF approach 
for the '"'Ca-l-'"'Ca reaction at the same energies as those in Fig.|4] 
Number of exchanged nucleon is superimposed by the solid cir- 
cles, solid triangles, and solid squares from high to low energies, 
respectively. 



energy of this system is 53.4 MeV. The magnitude of diffu- 
sion coefficient essentially depends on the size of the win- 
dow, the larger the window the larger the rate of change of 
nucleon exchange [i32J . At a given center-of-mass energy, 
the diffusion coefficient becomes maximum at the turning 
point where the size of the window is the largest. Also, be- 
cause of increasing overlap of the projectile and the target, 
the magnitude of the diffusion coefficient increases with 
energy. 

In Fig. |5] the variance of the fragment mass distribu- 
tions deduced from the SMF approach at the same energies 
as those in Fig. |4]are shown by lines. The corresponding 
evolution of the number of exchanged nucleons is superim- 
posed by the solid circles, solid triangles, and solid squares 
from high to low energies, respectively. The mass variance 
estimated from the SMF approach is consistent with this 
relation. We also estimate the variance of fragment mass 
distribution using the standard TDHF approach. In Table 1, 
the asymptotic values of the variances obtained from SMF 
and from TDHF are compared with those of the number of 
exchanged nucleons. The TDHF results are much smaller 
than the number of exchanged nucleons and are also much 
smaller than the results obtained from the SMF approach. 
The failure of the TDHF theory on the description of the 



Table 1. Asymptotic values of the variances obtained from TDHF 
(cr|pjjp) and SMF (cr\^) for the '*''Ca+'"'Ca reaction. Asymptotic 
values of the number of exchanged nucleons are also given in the 
last column. 
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variance of the fragment mass distribution has been recog- 
nized for a long time as a major limitation of the mean- 
field theory |11 12 13 1. It appears that the SMF approach 
cures this shortcoming of the mean-field theory. As seen 
from Fig.|5] not only the asymptotic value of cr^^ but also 
the entire time evolution is very close to the evolution of 

A^exc(f)- 



4 Nucleon drift in asymmetric system 

The work presented in the previous section has been ex- 
tended to asymmetric systems. Then, in addition to the dif- 
fusion coefficient, the drift term in Eq. ^ connected to the 
average nucleon transfer is different from zero. Its techni- 
cal details will be reported elsewhere ll35l . In order to ex- 
tend the SMF approach to the asymmetric system, includ- 
ing the motion of the neck position is necessary. Following 
the nucleon exchange picture |30|, the drift term can be 
estimated through the semiclassical expression of the SMF 
approach: 



VA(t) 



I 



dpx \Px - Pol 
2nti m 



(jco,/7.„0-/r(^o,/7x,0], (12) 



where po/w is the velocity of the neck position. On the 
other hand, the drift term in TDHF is given with the help 
of the Wigner function as 



VA{t) 



j 2nh 



■Po 



f{xo,px,t), (13) 



where the Wigner function is given from Eq. (O by 

fix, Px, = 2 ^P'^^^ + ^T^i^^ 0] . (14) 

The expression of Eq. ( fT2b can be seen as a semiclassical 
version of Eq. ( fT3l l. We have recently examined its validity 
by comparing the average number of nucleon transfer. 



Arit = oo) - Arit = 0) 



Jo 



dtVA(t), 



(15) 



where Arit) denotes the mass number of the target-like 
fragment at time t. 

Figure|6]shows the result of Eq. ( fTSl l obtained by TDHF 
and by the semiclassical expression of SMF as a function 
of center-of-mass energy for the head-on "^''Ca-n'^'^Zr reac- 
tion. The semiclassical expression which has been used at 
energies above the Coulomb barrier leads to values consis- 
tent with its fully quantal counterpart. However, it clearly 
breaks down at energies below the Coulomb barrier. This 
is due to the fact that the effect of nucleon tunneling is 
not treated well in the semiclassical expression below the 
Coulomb barrier, underlying the necessity of quantal de- 
scription of nucleon transfer 
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Fig. 6. The average number of nucleon transfer from '^''Ca to 
'"Zr evaluated from TDHF (solid line) and from the semiclassi- 
cal expression of SMF (dashed hne) are shown as a function of 
center-of-mass energy for the head-on "^^Ca+'^Zr reaction. The 
Coulomb barrier Vg is indicated as the arrow. 



5 Conclusion 

Mean-field dynamics and mean-field fluctuations using mi- 
croscopic time-dependent models are discussed in the con- 
text of low energy nuclear reactions. We have shown that 
the TDHF theory gives precise values of nucleus-nucleus 
potential and a universal behavior of energy dissipation. 
By projecting the SMF equation on the relative distance 
of colliding nuclei and the mass number of the projectile- 
Uke nucleus, we extract the corresponding diffusion coef- 
ficients. We have shown that the SMF approach correctly 
describes the mass variance of final fragments in transfer 
reactions at energies near the Coulomb barrier This gives 
a practical solution to properly describe mean-field fluctu- 
ations of one-body observables at low energies. 



This work is supported in part by US DOE Grant DE-FG05- 
89ER40530. 



References 

1. p. Bonche, S. E. Koonin, and J. W. Negele, Phys. Rev. 
C 13, (1976) 1226. 

2. S. E. Koonin, Prog. Part. Nucl. Phys. 4, (1980) 283. 

3. J. W. Negele, Rev. Mod. Phys. 54, (1982) 913. 

4. C. Simenel, D. Lacroix, and B. Avez, Quantum Many- 
Body Dynamics: Applications to Nuclear Reactions 
(VDM Verlag, Sarrebruck, Germany, 2010). 

5. K.-H. Kim, T. Otsuka, and P Bonche, J. Phys. G 23, 
(1997) 1267. 

6. C. Simenel, Ph. Chomaz, and G. de France, Phys. Rev. 
Lett. 86, (2001)2971. 

7. T. Nakatsukasa and K. Yabana, Phys. Rev. C 71, (2005) 
024301. 

8. J. A. Maruhn, P-G. Reinhard, P D. Stevenson, J. R. 
Stone, and M. R. Strayer, Phys. Rev. C 71, (2005) 
064328. 



EPJ Web of Conferences 



9. A. S. Umar and V. E. Oberacker, Phys. Rev. C 73, 
(2006) 054607. 

10. Y. Iwata, T. Otsuka, J. A. Maruhn, and N. Itagaki, 
Phys. Rev. Lett. 104, (2010) 252501. 

11. S. E. Koonin, K. T. R. Davies, B. Maruhn-Rezwani, H. 
Feldmeiner, S. J. Krieger, and J. W. Negele, Phys. Rev. C 
15, (1977) 1359. 

12. K. T. R. Davies, V. Marahn-Rezwani, S. E. Koonin, 
and J. W. Negele, Phys. Rev. Lett. 41, (1978) 632. 

13. C. H. Dasso, T. D0ssing, and H. C. Pauli, Z. Phys. A 
289, (1979) 395. 

14. Y. Abe, S. Ayik, P-G. Reinhard, and E. Suraud, Phys. 
Rep. 275,(1996) 49. 

15. D. Lacroix, S. Ayik, and Ph. Chomaz, Prog. Part. Nucl. 
Phys. 52, (2004) 497. 

16. R. Balian and M. Veneroni, Phys. Lett. B136, (1984) 
301; R. BaHan, P. Bonche, H. Flocard, andM. Veneroni, 
Nucl. Phys. A428, (1984) 79c. 

17. J. B. Marston and S. E. Koonin, Phys. Rev. Lett. 54, 
(1985) 1139. 

18. J. M. A. Broomfield and P. D. Stevenson, J. Phys. G 
35, (2008) 095102. 

19. C. Simenel, Phys. Rev. Lett. 106, (2011) 112502. 

20. S. Ayik, Phys. Lett. B658, (2008) 174. 

21. K. Washiyama and D. Lacroix, Phys. Rev. C 78, 
(2008) 024610. 

22. K. Washiyama, D. Lacroix, and S. Ayik, Phys. Rev. C 
79, (2009) 024609. 

23. S. Ayik, K. Washiyama, and D. Lacroix, Phys. Rev. C 

79, (2009) 054606. 

24. K. Washiyama, S. Ayik, and D. Lacroix, Phys. Rev. C 

80, (2009) 031602(R). 

25. J. O. Newton et al., Phys. Lett. B586, (2004) 219; 
Phys. Rev. C 70, (2004) 024605. 

26. V. Y. Denisov and W. Norenberg, Eur. Phys. J. A 15, 
(2002) 375. 

27. G. G. Adamian, R. V. Jolos, A. K. Nasirov, and A. 1. 
Muminov, Phys. Rev. C 56, (1997) 373. 

28. H. Esbensen, A. Winther, R. A. BrogUa, and C. H. 
Dasso, Phys. Rev. Lett. 41, (1978) 296. 

29. C. H. Dasso and R. Donangelo, Phys. Lett. B276, 
(1992) 1. 

30. H. Feldmeier, Rep. Prog. Phys. 50, (1987) 915. 

31. J. Randrup and W. J. Swiatecki, Ann. Phys. (N.Y.) 
125, (1980) 193; Nucl. Phys. A429, (1984) 105. 

32. J. Randrup, Nucl. Phys. A383, (1982) 468. 

33. H. Freiesleben and J. V. Kratz, Phys. Rep. 106, (1984) 
1. 

34. G. G. Adamian, A. K. Nasirov, N. V. Antonenko, and 
R. V. Jolos, Phys. Part. Nucl. 25, (1994) 583. 

35. B. Yilmaz, S. Ayik, D. Lacroix, and K. Washiyama, in 
preparation. 



